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Abstract 

We study charmless pure annihilation type radiative B decays within the QCD factorization approach. After adding 
the vertex corrections to the naive factorization approach, we find that the branching ratios of B® — > <j>y, B® p°y 
and Sj -¥ coy within the standard model are at the order of ^(1(T 12 ), ^(KT 10 ) and ^'(KT 11 ), respectively. The 
smallness of these decays in the standard model makes them sensitive probes of flavor physics beyond the standard 
model. To explore their physics potential, we have estimated the contribution of Z' boson in the decays. Within the 
allowed parameter space, the branching ratios of these decay modes can be enhanced remarkably in the non-universal 
Z' model: The branching ratios can reach to ^(10~ 8 ) for p°(ft))y and ^(1CT 10 ) for the B d — s> <j>y, which are 
large enough for LHC-b and/or Super B-factories to detect those channels in near future. Moreover, we also predict 
large CP asymmetries in suitable parameter space. The observation of these modes could in turn help us to constrain 
the Z' mass within the model. 
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1 Introduction 



Rare B decays induced by flavor changing neutral currents (FCNC) play important roles in particle physics, where they 
are always regarded as ideal places for probing signals of new physics. The GIM suppression of FCNC amplitude is 
absent in many new physics scenarios beyond the standard model (SM), which could give large enhancement of FCNC 
contributions over the SM predictions. However, due to our poor knowledge of non-perturbative QCD, predictions 
for many interesting exclusive decays are polluted by large hadronic uncertainties. Therefore, it would be of great 
interest to explore rare B decays, which are induced with few hadronic parameters as well as only by FCNC currents. 
Two body radiative B decays involve simple hadronic dynamics with only one hadron in the final states, so they suffer 
much less pollution than non-leptonic decays. 

In studying the radiative decays such as B — > K*y, p(o)y, the isospin breaking effects between the charged B and 
neutral B° in these modes are mainly from the annihilation type diagrams |[T] [2] [3j |4j . Many of pure annihilation type 
radiative decays, such as B — > <py and B — > J '/y/y, have been analyzed in the QCD factorization approach J5j|6l and in 
the perturbative QCD approach Q. We find the branching ratio of B — > 0y is at the order of 0(1O~ 11 ~ 10~ 12 ) in the 
SM. The decay rate is too small to be observed at presently running B factories, BaBar and Belle. Any measurements 
of the decays at BaBar and Belle would be direct evidences of new physics. In this work, we explore the decay 
B —> <py and similarly the decay B° s p°(o))y in the non-universal Z' model H), which could be naturally derived in 
certain string constructions (9|, E6 models ITOl and so on. Generally speaking, within the such model a flavor mixing 
can be induced at the tree level in the up-type and/or down-type quark sector after diagonalizing their mass matrices. 
In some new physics model, FCNCs due to Z' exchange can be induced by mixing among the SM quarks and the 
exotic quarks, which have been predicted to have different Z' quantum numbers. Here we will consider the model in 
which the interaction between the Z' boson and fermions are flavor non-universal for left handed couplings and flavor 
diagonal for right handed couplings. The effects of the Z' on other processes of the interest have been investigated in 
a number of papers such as lfTTl[T2ll . especially in B physics lfl~3l [141 [T5l [T6l . The recent review about Z' in detail is 
referred to Ref. ifTTl . 

To keep completeness, we first calculate these decays in the naive factorization approach. Then we add the vertex 
corrections to the four quark operators, which have been performed in the so called QCD factorization approach [18] 
in the SM, utilizing the light-cone wave functions of the light vector mesons. A similar work within the R-parity 
violating SUSY can be also found in Ref. 0. However, in this work we will revisit these processes with the updated 
parameters in the non-universal Z' model. 



2 Calculation in the Standard Model 



In the SM, the common starting point is the effective weak Hamiltonian which mediates flavor-changing neutral current 
transitions of the type b — > D (D = d,s) : 

G r / io.7y.8s 

•&ff = ~k L Wo ( ci o\ + c 2 o p 2 - v tb v; D £ qo 

V2 L p=u,c V ' z=3 
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Figure 1: The leading diagram for6 5 — > p°y. The symbol ® denotes the insertion of penguin operators O,. 



The explicit forms of the operators 0\ and the corresponding Wilson coefficients C, at the scale of jti = mj, can be found 

in Ref. |fl9l . V p (i)fc, V p ^ D are the Cabibbo-Kabayashi-Maskawa (CKM) matrix elements. According to the effective 

Hamiltonian ([T), we can draw the lowest order diagram of this channel, as shown in Fig. 1, which is dominated by the 

o 

photon radiated from light quark in B s meson. When the photon is radiated from heavy b quark, the energetic up or 
down quark will suppress the p° production by a power of A^ CD /m^, therefore, here we neglect its contribution. This 
point have been clearly discussed in B — ► K* y decays JUG]. 

For convenience, we denote that TJ* and £* are transverse polarization vectors of the final vector meson and 
photon, respectively. The photon energy and momentum are defined as Ey and q, and the momentum of B q meson is 
Pb = mv. In the rest frame of B q meson, we take the photon and the vector meson moving along the n_ = (1,0,0,-1) 
and n + = (1,0,0, 1) directions, respectively. Within the effective Hamiltonian and naive factorization hypothesis, we 
can write down the amplitudes as follows: 



A(B°^p<V) 



G F 



V^efpmpFy | -^vp^ferv-V + i [(t?1 ■ el)(v • q) - (tj! ■ q)(v ■ e* ± )} J , (2) 



A{B^ay) = ^- 



y^F e f p m p Fv | -e^p a ri*/e* ± v vPq a + i [{r\l ■ e* x ){v ■ q) - (r,* ± ■ q)(v ■ el)] | , (3) 



A(B° d -> 0y) = ^ 



- v tbV t * d I «3 + «5 - \a n - \-Cl 9 



x y^JpmpFv | -e^pa 777 el v vP q a + i [(ril ■ el)(v • q) - (ijl ■ q) (v ■ el)] I , (4) 



where a, is defined as the combination of the Wilson coefficients, 



at = Q 



Nr. 



(5) 
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Table 1 : Summary of input parameters 


A 


A P V A< ^ 4) t b° t b° A « a " 

MS s 


a s 


0.225 


0.818 0.141 0.348 250MeV 1.54ps 1.46ps 0.35 1/132 


0.214 


fs 


/b s fp fp f<j> fty fa 


fa 


216MeV 


236MeV 210MeV 150MeV 221MeV 175MeV 187MeV 


151MeV 




niB niB s m$ m p m m 






5.28GeV 5.36GeV 1.02GeV 0.77GeV 0.78GeV 





for an odd (even) value of i. The form factor Fy has been defined in Refs. f20l l2Tll22l 



{y(e*,q)\drn{l - 75)b\B° d (v)) = ^\%a e \ - F v (E r )e^ V p a e* v v p q a + iF A (E r )(e*q-v - q^ve*) 
In order to calculate the form factor Fy, we need two-particle light-cone projector for an initial B meson: 



(6) 



(7) 



where 4>£, (l + ) and 4>g 2 (Z + ) are the leading twist light-cone distribution functions |23j]. Thus, we obtain the standard 
result: 



F V (Ey) =F A (Ey 



Q.sf B M B r,, <MM 



dL 



U 



(8) 



dl + 

1+ ^B 

Consequently, we write down the helicity amplitudes for these channels as: 



2V2£ r 

where Q s = —1/3 is the charge of the s quark in units of the proton's charge. Because we have little knowledge about 
the distribution of the heavy meson, the integral in above formulae is often parameterized as : 

**(*♦■) 1 (9) 

(10) 

(ID 

(12) 
(13) 



B s ^py 



^B^ay 



'/++ 
l B d ^y 



i— ^Ana e F v f^m^M B 
i^f ^/\na e F y f^m^MB 



^B s ^py 



i^j= \J\na e F v f^m^MB 

•^B^ay — ^B d ^,$y 



3 3 
V ub V* s a 2 - V tb V* ( -fly + -ag 



2a 3 +2a 5 + -C 7 + ^a 9 



1 1 



v tbV t d fl 3 + as - ~a 7 - -a g 



0. 



Depending on the parameter values listed in Table 1, one can get the averaged branching ratios as: 



^(B?->p°y) = 1.1 x lO" 10 ; 




3§{B°^ coy) = 5.6 x 10~ n ; 




^(B^</>°y) = 1.7xl0- 13 . 


(14) 



Within the naive factorization hypothesis, because of no strong phases entering into these processes, there should not 
exist any CP asymmetry for the processes. 




(a) 




(b) 





Figure 2: Non-factorizable diagrams for/J^ — > p°y. 



Up to now in our calculation, non-factorizable contributions have been neglected. As next step, we add the vertex 
corrections and the leading non-factorizable diagrams, shown in Fig. 2. To achieve the goal, the QCD factorization 
framework |fl8ll proposed by Beneke, Buchalla, Neubert and Sachrajda is very suitable to be applied. To calculate 
non-factorizable diagrams, we also need the two-particle light-cone projector of the vector mesons: 



fp m P 



(15) 



where g^(u) and g^"\u) are twist-3 distribution amplitudes of vector mesons, and explicit formulae can be found in 
Ref. El. 

After adding the contributions, the form of amplitudes for the decay modes becomes the similar, just replacing a, 
by a-, which involve the ff(a s ) corrections, a-'s are calculated to be 

, a s C F fp 

° 2 = ai + TnT c T p ClFu 

, a s C F fp ni7 

° 3 = ai + VnW c T p C ^ 

a s C F fp 
fl5 = a5 + ^T c Tp 6 2 ' 

i . a s C F fp 

a 7 = a 7 + — ——C%F 2 , 
AnN c fp 

. u s C F fp 



Clg 



(16) 
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Table 2: The values of a, in different scenario. In the LO column, a, are defined in Eq. ©; in the NLO column, the 
values are a\ defined in Eq. (fT6] i; in the Z' Model, they are values of a' t + Aa'j, which are defined in Eq. ( f29l >. 





LO 


NLO 


Z' Model 


a 2 


0.170 


0.149-/0.010 




0.149- 


-iO.010 




a 3 


0.002 


-0.002 - zO.002 


(0.019 4 


-(0.007)^ 


+ (-0.002 


- /0.002) 


as 


-0.005 


0.003 + j'0.002 


(0.009 


-(0.008)£e' 


> + (0.003 -+ 


(0.002) 


a-, 


0.000 


-0.000 + /0.000 


(3.798- 


-/0.067)4«''* 


+ (-0.000 


fiO.OOO) 


ag 


-0.008 


-0.008 - /0.000 


(3.932- 


(0.050)<^ 


+ (-0.008 


- /0.000) 



where Fi^ arise from one gluon exchange between the two currents of color-octet operators as shown in Fig. 2, 

J" 



Fi = 



, / gj (") (v), x 

„ rfM 1 - 4 *x W 



-14-3(7T-121n — 



ilnM- — +2Li 2 ( ) 

1 — M / 3 M 



(17) 



F 2 



£du (g ( l\u)- 



sY (") 



14-3/7r-121n 



u \ it „ . ,u — 1 . 
InM- — +2Li 2 ( 

1 — U J 3 M 



(18) 



Here we have neglected the small effect of box diagrams and the diagrams with photon radiating from energetic light 
quarks, which are further suppressed by Aqcd/Mb- Including &{a s ) contributions, the averaged branching ratios in 
the SM are estimated to be 



&(B° S -> coy) 



1.1 x 10 



-10 



2.3 x 10 



^(B° d -> <j>y) = 2.9 x 10 



-12 



(19) 



Comparing with the results in Eq. (TBI of the naive factorization, one finds the branching ratio of B® — > p°y almost 
unchanged. To find out the reason why the correction dose not take an effect, we list the values of a\ of these decay 
modes in the Table 2. From the table, we find that the corrections to a-i .9 are very small and can be neglected. 
Although Wilson coefficients of QCD penguin operators changed a little, but they give no contribution because the 
quark component of p° is (uu — dd)/y/2. For a 2 , it changes much, but the correction is suppressed by the CKM 
elements. So, the unchanged branching ratio is quite reasonable. As for B® — > coy, the decrease of as can cause that 
the branching ratio becomes even smaller than that of the naive factorization. For the decay B® — > (j>y, the increase of 
the ratio mainly comes from the change of a 3 and as. 

Because there are both weak and strong phases in the decay modes B® — > p°y and B® — > coy, we can get the CP 
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asymmetries of these two channels as follows, 



?->py) = 3%; 
° -> coy) = -27%, 



(20) 



by the definition of CP asymmetry 



BR(B° -> Vy) — BR(B® -> Vy) 
BR(B° -> Vy) +BR(B° -> Vy) ' 



(21) 



For the decay mode — ► 07> there is only weak phase from V t bV* d , so that the CP asymmetry in this decay disappears 
within the SM. 



3 Calculation in the Non-universal Z f Model 



Now we consider the effects due to an extra 1/(1)' gauge boson Z'. Usually, the flavor mixing can be induced at the 
tree level in up-type and/or down-type quark sector after diagonalizing their mass matrices. In some new physics 
model, FCNCs due to Z' exchange can be induced by mixing among the SM quarks and the exotic quarks, which is 
predicted and can make different Z' quantum numbers after the mixing. Here we will consider the model in which the 
interaction between the Z' boson and fermions are flavor non-universal for left handed couplings and flavor diagonal 
for right handed couplings. For simplicity, we neglected the mixing between the Z° and Z' and the evolution effect 
from the high scale M z > to the M\y scale. 

We start to set up the relevant interactions with the new Z' gauge particle. Following the convention in Ref. J8), 
we write the couplings of the Z' -boson to fermions as 



(22) 



where i is the family index and y labels the fermions and Pl.r = (1 =F Ys) /2. According to some string construction 
or GUT models such as £g, it is possible to have family non-universal Z' couplings. That is, even though E,' R are 
diagonal, the couplings are not family universal. After rotating to the physical basis, FCNCs generally appear at tree 
level in both left handed and right handed sectors. Explicitly, 



B^ = V VL e^v; L , B*« =V VR e**v; R . 



(23) 



Moreover, these couplings may contain CP-violating phases beyond that of the SM. The effective Hamiltonians de- 
scribing the transition mediated by the Z' boson have the form as: 



V/ ■ 



e f f (b^sqq) = -^V th V* 



Jt?f f (b -> dqq) 



V2 

4G F 
"72 



V,bV*i 



g'M z 
giM z ,J 

2 



' (B L m 9 +B«0 7 ) 



( g'M z \ 



\g\M Z ' J V tb V t 



» (B L o 9 +B R n 7 ) 



td 



(24) 
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L(R) 

where g\= e/ (sin 9w cos 9w) and Bjj denote the left (right) handed effective Z' couplings of the quarks i and j at 
the weak scale. The diagonal elements are real due to the hermiticity of the effective Hamiltonian but the off diagonal 
elements may contain effective weak phases. With the definition 

^(^) 2 , (25) 
giM z i 

we can parameterize these coefficients as 



V V,bV td J V ' 6 ' d J 

where (j)' = <p — /3, (0 is the weak phase associated with B^ lb ). 

l(r\ L( R} l,( R) 

In the following discussion, we adopt BT ' * -2B2;> * -2BT> for convenience (note that a possible negative 



sign can be accounted for by shifting <j) by n), which has been stated in detail in Ref. ITOI . In order to see the effect 
of Z' boson, we have to know the values of the AC7 and AC9 or equivalently B^j b , B^ b and Bqf. Generally, one expects 
g'/gi ~ 1, if both the U(l) gauge groups have the same origin from some grand unified theories, and Mz/M z i ~ 0.1 
for a TeV scale neutral Z' boson, which yields y ~ 10~ 2 . However, in Ref. |fl3ll assuming a small mixing between 
Z — Z' bosons the value of y is taken as y ~ 10~ 3 . It has been shown in lfT4l that the mass difference of B s —B s 
mixing can be explained if \Bj b \ ~ |^V^*|. Similarly, the CP asymmetry anomaly in B — > 7r7T can be resolved if 
\B L sb B L s f\ ~ |V,iV r *|. So, we assume that 

111 I = 1^1 = llfl = \\tt U \ 6 (10- 3 ,10- 2 ). (27) 

Assuming only left handed couplings are present, the bound on FCNC Z' coupling (B^ lb ) fromB — B° mass difference 
has been obtained in fl4l as 

y|Re(B^) 2 |<5xl0- 8 , y|Im(z4) 2 |<5xl(T 8 . (28) 

Using y ~ 10~ 2 , one can obtain a more stringent bound on \Bj b \ < 10~ 3 . From these two relations one can obtain 
\B^ S \ ~ 1. Thus, it is expected that % 2 ~ 10~ 3 with the CKM matrix elements considered. However, in our analysis 
here we vary their values within the range \^\ & (10~ 3 , 10~ 2 ), since the major purpose of this work is searching for 
new physics signal rather than obtaining acute numerical results. 

It is noted that the other Wilson coefficients may also receive contributions from the Z' boson through renormal- 
ization group (RG) evolution. With our assumption that no significant RG running effect between M' z and Mw scales, 
the RG evolution of the modified Wilson coefficients is exactly the same as the ones in the SM lfl9l . Using the values 
of these coefficients at nij, scale we can analogously obtain the new contribution to the transition amplitude as done in 
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10 

Figure 3: For the decay mode B — > <j>y, variation of the CP averaged branching ratio (in units of 10~ 10 ) with E, (in units 
of 10~ 3 ) and the new weak phase (left panel) and the variation of direct CP asymmetry (in %) with the new weak 
phase (right panel) where the solid, dot-dashed and dashed lines correspond to E, = 0.001 , 0.005 and 0.01 . 
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Figure 4: For the decay mode B s —> p°j, variation of the CP averaged branching ratio (in units of 10~ 8 ) with t, (in 
units of 10~ 3 ) and the new weak phase (left panel) and the variation of direct CP asymmetry (in %) with the new 
weak phase <p (right panel) where the solid, dot-dashed and dashed lines correspond to £, = 0.001 , 0.005 and 0.01 . 
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30 



( = 0.001 




Figure 5: For the decay mode B s — > oy, variation of the CP averaged branching ratio (in units of 10~ 8 ) with B, (in 
units of 10~ 3 ) and the new weak phase (left panel) and the variation of direct CP asymmetry (in %) with the new 
weak phase (right panel) where the solid, dot-dashed and dashed lines correspond to ^ = 0.001 , 0.005 and 0.01 . 



the case of Z boson. The Aa^ induced by the Z'are given as: 



Aa' 2 
Aa' 3 



Aa' 7 
Aa' 9 



ACi 
AC 3 
AC 5 
AC 7 - 
AC 9 



AC 2 
AC 4 



Afl', = AC 5 + — - 



AQ 
3 



AQ 
3 

ACio 



«v C F fp 

^ f -tAC 4 F h 
cc s C F fp „ 

AnNcfp 
AnNcfp 



AC l0 F u 



and the contributions of new physics can be formulated as 





G F 


^t% 7 = 




A ^io r = 


.G F 
l ~ 




.G F 
l V=2 



s/Ana e F v f^m^MB 



V uh V: s Aa' 2 -V tb V; s 
V ub V: s Aa' 2 -V tb V* 



-Ml + ~ Aflg 



2Aa' 3 + 2A«5 + - Aa'-j + -Aa' g 



1 



1 



-V tb V* d ( Aa' 3 + Aa' 5 - -Aa' 7 - -Aa' 9 



(29) 

(30) 
(31) 
(32) 



Noted that Aajs involve the new weak phase, which may change the CP asymmetries remarkably. Now using |^i | = 
|^2| = % and taking the decay mode B s —> p°y as an example, we also list the correction to d i from Z' boson in the 
last column of Table 2. From the table, we note that the corresponding Wilson coefficient of electro-weak penguin 
is enhanced remarkably with suitable parameter which may affect the branching ratio and other observed values. 
In Fig. [3] Fig. [4] and Fig- H we show the variation of the CP averaged branching ratios of B — s- (j)y, B s —> p°y and 
B s coy with £, and the new weak phase 0" (left panel) and the corresponding direct CP violation with 0" (right 
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panel), respectively. As anticipated, if the E, = 0.01, the branching ratios can be enhanced remarkably, which can 
reach to <^(10~ 8 ) forB s — > p («)y and ^?(10~ 10 ) for the Bj — > (j)°y. All results are enhanced two orders of magnitude 
over the predictions of the SM. From the figures, we can also argue that there have been significant enhancements in 
the branching ratios for large ^ , or in other words for a lighter Z' boson. In the experimental side, these results may 
be inaccessible at the Belle and BaBar presently. However, it is large enough for LHC-b and/or Super B-factories. 
Moreover, we find that these decay modes may have large CP asymmetries when % = 0.001 and suitable weak phase (j>. 
It implies that the contributions from new physics and from the SM can be comparable, and the interference between 
them leads to large asymmetries. Furthermore, future observations of these modes could in turn help us to constrain 
the mass of Z' boson within the model. 



4 Conclusion 

In this work, we have studied pure annihilation type radiative processes lf d — > <j>y, B s — > py and — > coy within 



the QCD factorization. After adding the vertex corrections to the naive factorization approach, we find that the non- 
factorizable contributions can enhance the branching ratio of B® — >• y decays, however, the branching ratio of B® —> p y 
is almost unchanged, but for B® — > coy the branching ratio is even lowered, because the corrections to EW penguin 
operators are much smaller than those to QCD penguin operators. The smallness of these decays within the SM makes 
them sensitive probes of flavor physics beyond the SM. To explore new physics potential, we have estimated the 
contribution of the non-universal Z' model to the decays. If £ = 0.01, the branching ratios can be enhanced remarkably, 
and reach to <^(10~ 8 ) for B s -» p°(co)y and ^(10 -10 ) for the B c j — » <j)°y. Moreover, we have also predicted large CP 
asymmetries in suitable parameter spaces. These results can be tested at the LHC-b and/or Super B-factories in future. 
The observations of these modes could in turn help us to constrain the mass of Z' within the model. 

Acknowledgments 

The work of C.S.K. was supported in part by Basic Science Research Program through the NRF of Korea funded by 
MOEST (2009-0088395) and in part by KOSEF through the Joint Research Program (F01-2009- 000-10031-0). The 
work of Y.Li was supported by the Brain Korea 21 Project and by the National Science Foundation under contract 
Nos. 10805037 and 10625525. Y.Li thanks Prof Y.-D Yang for valuable discussion. 

References 

[1] A. Ali and A. Parkhomenko, |arXiv:hep-ph/06 1 1 49] 



[2] S. W. Bosch and G. Buchalla, Nucl. Phys. B 621, 459 (2002) | |arXiv:hep-ph70T0 6081|. 

[3] P. Ball, G. W. Jones and R. Zwicky, Phys. Rev. D 75, 054004 (2007) |arXiv:hep-ph/06 12081). 



11 



[4] W. Wang, R. H. Li and C. D. Lu JarXiv:07lT0432 [hep-ph]. 

[5] X. Q. Li, G. R. Lu, R. M. Wang and Y. D. Yang, Eur. Phys. J. C 36, 97 (2004) | |arXiv:hep-ph/03052831 . 
[6] G. R. Lu, R. M. Wang and Y. D. Yang, Eur. Phys. J. C 34, 291 (2004) | |arXiv:hep-ph/0308256) . 
[7] Y. Li and C. D. Lu, Phys. Rev. D 74, 097502 (2006) | |arXiv:hep-ph/0605220) . 
[8] P. Langacker and M. Plumacher, Phys. Rev. D 62, 013006 (2000) larXiv:hep-ph/0001204) . 
[9] G. Buchalla, G. Burdman, C. T. Hill and D. Kominis, Phys. Rev. D 53, 5185 (1996) | |arXiv:hep-ph/95 10376] . 
[10] E. Nardi, Phys. Rev. D 48, 1240 (1993) larXiv:hep-ph/9209223| . 

[11] J. Erler, P. Langacker, S. Munir and E. R. Pena, JHEP 0908, 017 (2009) larXiv:0906.2435i rhep-phll. 
[12] S. L. Chen and N. Okada, Phys. Lett. B 669, 34 (2008) HarXiv:0808.0331i rhep-ph11. 

[13] V. Barger, C. W. Chiang, P. Langacker and H. S. Lee, Phys. Lett. B 598, 218 (2004) |arXiv:hep-ph/0406126| ; 
V. Barger, L. Everett, J. Jiang, P. Langacker, T. Liu and C. Wagner, Phys. Rev. D 80, 055008 (2009) 
llarXiv:0902.4507l [hep-ph]]; V. Barger, L. L. Everett, J. Jiang, P. Langacker, T. Liu and C. E. M. Wagner, 
larXiv:0906.3745l rhep-ph] : V. Barger, C. W. Chiang, P. Langacker and H. S. Lee, Phys. Lett. B 580, 186 (2004) 
QarXiv:hep-ph/0310073| . 

[14] K. Cheung, C. W. Chiang, N. G. Deshpande and J. Jiang, Phys. Lett. B 652, 285 (2007) |arXiv:hep-ph/0604223) . 
C. W. Chiang, N. G. Deshpande and J. Jiang, JHEP 0608, 075 (2006) | |arXiv:hep-ph/0606122) . C. H. Chen and 
H. Hatanaka, Phys. Rev. D 73, 075003 (2006) | |arXiv:hep-ph/0602140) . 

[15] Q. Chang, X. Q. Li and Y. D. Yang, JHEP 0905, 056 (2009) flarXiv:0903.0275i rhep-phll; Q. Chang, X. Q. Li and 
Y. D. Yang. larXiv:0907.4408l rhep-ph] . 

[16] C. H. Chen, |arXiv:091 1.34791 [hep-ph] . C. W. Chiang, R. H. Li and C. D. Lu, larXiv:091 1.23991 [hep-ph] . R. Mo- 
hanta and A. K. Giri, Phys. Rev. D 79, 057902 (2009) narXiv:08 12. 1 8421 [hep-ph] 1 . 

[17] P. Langacker. larXiv:0801 . 1 3451 fhep-phl . 

[18] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999) 
QarXiv:hep-ph/9905312| ; M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nucl. Phys. B 591, 313 
(2000) | |arXiv:hep-ph/0006124) . 

[19] For a review, see G. Buchalla, A.J. Buras, M.E. Lautenbacher, Rev. Mod. Phys. 68, 1 125 (1996). 

[20] E. Lunghi, D. Pirjol and D. Wyler, Nucl. Phys. B 649, 349 (2003) | |arXiv:hep-ph/0210091) . 

[21] S. Descotes-Genon and C. T. Sachrajda, Nucl. Phys. B 650, 356 (2003) flarXiv:hep-ph/0209216) . 

[22] P. Ball and E. Kou, JHEP 0304, 029 (2003) | |arXiv:hep-ph/0301135] . 



12 



[23] A. G. Grozin and M. Neubert, Phys. Rev. D 55, 272 (1997) | |arXiv:he p-ph/9607366|. 



[24] P. Ball, V. M. Braun, Y. Koike and K. Tanaka, Nucl. Phys. B 529, 323 (1998) [arXiv:hep-ph/9802299|. 



13 



